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ABSTRACT 

We present nuclear spectral energy distributions (SEDs) in the range 0.4 — 16 /im for an expanded CfA 
sample of Seyfert galaxies. The spectral indices (/„ oc i/ _QIR ) from 1 — 16 /im range from am ~ 0.9 to 3.8. 
The shapes of the spectra are correlated with Seyfert type in the sense that steeper nuclear SEDs (vf v 
increasing with increasing wavelength) tend to be found in Seyfert 2s and flatter SEDs (yf v ~ constant) 
in Seyfert 1 — 1.5s. The galaxies optically classified as Seyferts 1.8s and 1.9s display values of am as 
in type 1 objects, or values intermediate between those of Seyfert Is and Seyfert 2s. The intermediate 
SEDs of many Seyfert 1.8 — 1.9s may be consistent with the presence of a pure Seyfert 1 viewed through 
a moderate amount (Ay <J 5 mag) of foreground galaxy extinction. We find, however, that between 10 
and 20% of galaxies with broad optical line components have steep infrared SEDs. 

Torus models usually adopt high equatorial opacities to reproduce the infrared properties of Seyfert 
Is and 2s, resulting in a dichotomy of infrared SEDs (flat for type Is, and steep for type 2s). Such a 
dichotomy, however, is not observed in our sample. The wide range of spectral indices observed in the 
type 2 objects, the lack of extremely steep SEDs, and the large numbers of objects with intermediate 
spectral indices cannot be reconciled with predictions from existing optically thick torus models. We 
discuss possible modifications to improve torus models, including low optical depth tori, clumpy dusty 
tori, and high-optical-depth tori with an extended optically thin component. 

Subject headings: Galaxies: active — Galaxies: nuclei — Galaxies: photometry — galaxies: Seyfert — 
infrared: galaxies 



1. INTRODUCTION 

The discovery that Seyfert 2 galaxies such as NGC 1068 
can have reflected or polarized broad-line emission has led 
to a 'unified' model in which apparent differences in the 
properties of active galactic nuclei (AGNs) are interpreted 
in terms of viewing angle toward a similar intrinsic un- 
derlying structure (Antonucci 1993). In this unification 
paradigm, an optically thick circumnuclear dusty torus ab- 
sorbs a significant fraction of the optical/UV/X-ray lumi- 
nosity of the active nucleus and reradiates this energy at 
infrared wavelengths. One of the challenges of the unified 
model has been to reconcile the observed infrared emission 
with that predicted from the absorbing torus. 

Important observational constraints for the torus models 
come from the shape of the infrared spectral energy distri- 
butions (SEDs) and the silicate feature at 9.7 /im. Previ- 



ous studies have shown that Seyfert galaxies are stronger 
radiators in the mid-infrared than most non-Seyfert galax- 
ies (e.g., Spinoglio et al. 1995; Krabbe, Boker, & Maiolino 
2001) and that their infrared continua have a range of 
shapes from nearly flat {vf v ~ constant) to steeply rais- 
ing with increasing wavelengths in the region between 1 
and 10 /im. Seyfert 2 galaxies generally have SEDs steeper 
than Seyfert 1 galaxies (e.g., Rieke 1978; Edelson, Malkan, 
& Rieke 1987; Ward et al. 1987; Fadda et al. 1998; Alonso- 
Herrero et al. 2001). Torus models predict SEDs that are 
strongly dependent on the orientation of the torus with 
respect to the line of sight of the observer, as well as on 
the density distribution and geometry of the dusty absorb- 
ing material (Pier & Krolik 1993; Granato & Danese 1994; 
Efstathiou & Rowan-Robinson 1995; Granato, Danese, & 
Franceschini 1997). The emission between 1 and 10 /im 
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is thought to arise from hot dust (up to ~ 1000 K) (e.g., 
Rieke & Lebofsky 1981; Barvainis 1987) in the inner part 
of the torus and hence hidden from sight in edge-on sys- 
tems that appear to be of type 2. To match the 9.7 /im 
silicate feature, apparently often in absorption in Seyfert 
2 and not present or slightly in emission in Seyfert 1 spec- 
tra (Rieke & Low 1975a; Roche et al. 1991; Clavel et al. 
2000), the tori are assumed to be optically thick even to 
far- infrared wavelengths. 

Studies of AGNs have often focused on high luminos- 
ity objects (QSOs and bright Seyfert Is) since in these 
objects the active nucleus dominates the emission of the 
host galaxy. These studies avoid the complications of in- 
terpretation when the host galaxy emission is a signifi- 
cant fraction of the total (e.g., Kotilainen et al. 1992; 
Alonso-Hcrrcro, Ward, & Kotilainen 1996). Because mod- 
elling of SEDs is only available for a handful of relatively 
high- luminosity nearby AGNs (e.g., Efstathiou, Hough, & 
Young 1995; Alexander et al. 1999; Ruiz et al. 2001), the 
torus parameters arc not well constrained. The high angu- 
lar resolution of the Hubble Space Telescope (HST) now 
allows us to probe the nuclei with a beam area about 30 
times smaller than is typically achieved with ground-based 
observations at visible and near-infrared wavelengths. Im- 
agers in the mid-infrared also provide a significant im- 
provement in angular resolution over previous studies. 
Such data enable us to separate the nuclear emission from 
that of the surrounding galaxy with unprecedented accu- 
racy as well as to probe for lower luminosity active nu- 
clei. Recent imaging from the Infrared Space Observatory 
(ISO) also allows us to measure the nuclear emission in 
the mid-infrared with improved sensitivity compared with 
previous observations. 

We have used these capabilities on an expanded ver- 
sion of the CfA sample of Seyfert galaxies. The infrared 
SEDs of the CfA sample have been studied by Edelson et 
al. (1987), Pier & Krolik (1993), Maiolino et al. (1995) 
and Perez Garcia & Rodriguez Espinosa (2001). Fadda 
et al. (1998) compiled data from the literature and IRAS 
to produce infrared SEDs of a heterogeneous sample of 
Seyfert galaxies, but they mainly focused their efforts on 
fitting infrared colors to their model outputs. All these 
works noted that the relatively large apertures used were 
likely to be contaminated by emission from the host galax- 
ies, a particular problem for the Seyfert 2 galaxies. To 
address this concern, in this paper we present high reso- 
lution IRTF K' (2.1 /im) and L (3.5 /im) band imaging 
of the CfA Seyfert galaxies, from which we have esti- 
mated the non-stellar nuclear fluxes. We combine these 
fluxes with measurements of the unresolved nuclear emis- 
sion from optical HST/WFPCl and HST/WFPC2 and 
near-infrared iKT/NICMOS and mid-infrared measure- 
ments made from ISOCAM images, and/or small-aperture 
ground-based 10 ^im measurements from the literature. 
This suite of data is used to compare the properties of 
the non-stellar SEDs with the predictions of models for 
dusty tori. 

2. THE SAMPLE 

2.1. The CfA sample 

Tables 1 and 2 (see also Section 2.2) list the sample of 
galaxies on which this paper is based. The galaxies are 



from the CfA redshift survey, which is drawn from the 
fraction of the sky defined either by <5 > 0° and b > 40° 
or (5 > -2°. 5 and b < -30°, and m Zw < 14.5. The 
AGNs were identified by means of optical spectroscopy, 
thus defining a "complete" magnitude-limited sample of 
49 Seyfert galaxies (Huchra & Burg 1992). Because it 
is not color-selected, this sample should be relatively free 
of selection effects that tend to enhance the proportion 
of galaxies with anomalously strong emission in the color 
used for selection (Huchra & Burg 1992; Osterbrock & 
Martel 1993). 

Spectroscopic classifications for the CfA survey are 
taken from Osterbrock & Martel (1993) and references 
therein. However Ho, Filippcnko, & Sargent (1997) iden- 
tified 5 Seyfert 2 galaxies as having lower Seyfert types 
than listed by Osterbrock & Martel (1993). In these cases 
we list in Table 1, column 2 the identification from Os- 
terbrock & Martel (1993) on the left and that from Ho et 
al. (1997) on the right. Ruiz, Rieke, & Schmidt (1994) de- 
tected broad components in near-infrared lines in Mkn 334 
and NGC 5252 (He i at 1.083 /im) and Mkn 533/NGC 7674 
(He I and Pa/3). Optical polarized broad-line components 
were detected by Young et al. (1996) in NGC 4388, 
NGC 5252 and Mkn 533. Tran (2001) has found a hid- 
den broad line region in NGC 7682. For the Seyfert 2s in 
Table 1 we also indicate whether such hidden broad-line 
regions have been detected. 

2.2. The extended CfA sample of Seyfert galaxies and its 

completeness 

The CfA sample is usually regarded as a complete 
sample of optically selected Seyfert galaxies. However 
Maiolino & Rieke (1995) show that it is incomplete both 
in terms of low luminosity nuclei and edge-on galaxies. Ho 
& Ulvestad (2001) also discuss selection biases that may 
affect the completeness of CfA sample. Because of the 
method used to identify AGNs in the CfA sample, only 
Seyfert 2s having the broadest, most extended wings in the 
observed emission lines were included in the final sample. 
The distinction in the original classifications between LIN- 
ERs and Scyfcrts is also not rigorous. According to these 
authors the CfA sample is likely to be complete only for 
relatively bright objects, and the Seyfert 2s in particular 
are intrinsically more luminous and more radio-powerful 
than is typical. 

To alleviate in part the possible biases in the original 
CfA sample, we have defined an extended sample. It in- 
cludes all the Seyfert galaxies in the original CfA sample 
as well as those galaxies initially classified as LINERs by 
Huchra & Burg (1992), but later reclassified as Seyfert 
galaxies by Ho et al. (1997) from high S/N optical spectra. 
We therefore add a total of nine galaxies (see Table 2) to 
the original CfA sample, thus increasing its size by ~ 20% 
(a total of 58 galaxies). 

Our sample is still probably weighted against low- 
luminosity type 2 Seyferts. Table 3 compares the fractions 
of type 1 (1 — 1.5) and type 2 (1.8 — 2), as classified from 
optical spectroscopy in a number of nearby AGN samples: 
1.) the RSA sample (see Maiolino & Rieke 1995); 2.) the 
Palomar survey (Ho et al. 1997); and 3.) the infrared- 
selected 12 fxm sample (Spinoglio & Malkan 1989). Since 
the CfA sample shows a deficit of highly inclined galaxies 
(see Maiolino & Rieke 1995, and Section 5.5), in this table 
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we make the same comparison for galaxies with axial ra- 
tios b/a > 0.6. Galaxy inclinations are taken from Ho et 
al. (1997) and from the RC3 (de Vaucouleurs et al. 1991). 
We find that in low inclination galaxies, both the extended 
CfA sample and the 12 /im sample have slightly smaller 
fractions of type 2s than the Palomar and RSA samples, 
although the numbers are not different at a statistically 
significant level. If, as claimed by Ho & Ulvestad (2001), 
the CfA sample only contains Seyfert 2s with relatively 
bright broad emission lines, then the CfA sample would 
be biased against low luminosity type 2 AGNs rather than 
against pure type 2 objects (that is, those with an equa- 
torial or near equatorial view of the AGN). And indeed, 
it appears from Table 3 that the Palomar survey has the 
highest fraction of type 2 objects, probably because the 
high quality of its spectra allows detection of very low lu- 
minosity objects. 

3. OBSERVATIONS 
3.1. IRTF K' and L band observations 

We have obtained high resolution K' (A ccn trai = 
2.12 /im) and L band (A con t r ai = 3.51 /im) imaging of 34 
CfA Seyfert galaxies using NSFCam (Rayner et al. 1993; 
Shure et al. 1994) at the 3.8 m NASA IRTF telescope on 
Mauna Kea. The observations were made in August 1999, 
March 2000, August 2000 and May 2001. We used the 
0.3"/pixel plate scale. The nuclei of the galaxies were suf- 
ficiently bright that we were able to use the tip-tilt capabil- 
ity during the first three campaigns (it was not available 
during the fourth campaign), thus improving the image 
quality. Typical integration times were 300 seconds in the 
K' band and 720— 1200 seconds (depending on the galaxy 
brightness) in the L band. In Figure 1 we show three 
galaxies displaying extended emission in the L band. The 
spatial resolution (FWHM) in the L band as measured 
from the standard stars was 0.6" — 0.9". 

Standard data reduction procedures were applied for the 
K' images and for bright sources in the L band. Because 
of the short individual exposures in the L band (to avoid 
saturation) , a number of galaxies were not detected on in- 
dividual exposures, and thus the images had to be shifted 
to a common position and coadded using the nominal tele- 
scope offsets (blind offsets, see Table 4). In these cases the 
morphology and flux of the unresolved component cannot 
be determined, and only the integrated flux can be mea- 
sured. 

Conditions were photometric during all runs except for 
one night during the August 2000 campaign. Through- 
out the nights we obtained observations of standard stars 
from Elias et al. (1982) for photometric calibration. The 
typical rms scatter in these calibration observations is 
0.04 — 0.06 mag in both K' and L bands. The major source 
of uncertainty is the background subtraction, especially at 
the longer wavelengths (L) for the faint sources. In Table 4 
we give the K' and L band aperture photometry. 

3.2. Optical and near-infrared wavelengths: HST archival 

images 

iTST/NICMOS images of the CfA sample were taken 
for 22 Seyfert 1.8 — 2 galaxies with the NIC1 cam- 
era (0.043"/pixel) at 1.1 /zm (F110W filter) and 1.6 /im 
(F160W filter) as part of proposal 7867 and are de- 



scribed by Martini & Pogge (1999) and Martini et al. 
(2001). NGC 1068 was a guaranteed time target. In 
addition, some of the CfA Seyfert 1 — 1.5 galaxies were 
observed as part of a number of GTO and GO propos- 
als, mainly through the NIC2 (0.076"/pixel) F160W fil- 
ter. The F110W and F160W filters are roughly equivalent 
to the ground-based broadband J and H filters. Details 
on the ffST/NICMOS observations, data reduction and 
fluxes for the CfA Seyferts can be found in Quillen et al. 
(2001). 

When possible we also used WFPC2 archival images 
through the F606W filter (described by Malkan, Gorjian, 
& Tarn 1998). If these images were saturated we used im- 
ages taken through other filters or obtained the fluxes for 
the unresolved component from the literature (Ho & Peng 
2001 for WFPC1 or WFPC2). See Tables 1 and 2 for the 
filters. 

3.3. Mid-infrared wavelengths 

Mid-infrared (6.75 and 9.63 /im) fluxes measured from 
/SO/ISOCAM images were taken from Clavcl et al. 
(2000). The FWHM of the ISOCAM PSF is approximately 
^4 — 5" (Clavel et al. 2000). For a few galaxies which 
were not part of this program we obtained ISOCAM re- 
duced images from the ISO archive. We also make use 
of the 16 /im /SO/ISOPHOT fluxes measured for the CfA 
Seyferts by Perez Garcia & Rodriguez Espinosa (2001). 
In addition, we obtained ground-based small aperture N- 
band (10.6 /im) fluxes from the literature, both to compare 
with the lower resolution ISOCAM observations, and to 
use for those galaxies where ISOCAM data are not avail- 
able. 

In Table 1 and Table 2 we summarize the available 
iKT/WFPCl, WFPC2 and NICMOS, IRTF K'L and ISO 
and ground-based mid-infrared observations for the CfA 
Seyfert galaxies, along with appropriate references when 
the data were taken from the literature. 

4. THE UNRESOLVED EMISSION 

The contribution from the underlying galaxy can dom- 
inate the central emission in Seyfert 2 galaxies, especially 
at wavelengths shorter than approximately 2 /im (see e.g., 
Alonso-Herrero, Ward, & Kotilaincn 1996; Alonso-Herrero 
et al. 2001) and can also be significant in Seyfert 1 galax- 
ies (e.g., Kotilainen et al. 1992). It is therefore essential 
to eliminate the stellar contributions to isolate the non- 
stellar emission from the AGN. Since the data used in this 
work have been obtained with different instruments, the 
methods for obtaining the unresolved emission are slightly 
different, as we describe below. 

4.1. HST images 

In Quillen et al. (2001) we estimated the flux from the 
unresolved component in the F160W filter (1.6 /im) by fit- 
ting the sum of a galaxy profile and the point spread func- 
tion (PSF) for a large number of Seyfert galaxies. For the 
underlying galaxy the profiles were fitted by exponential 
and power laws convolved with the appropriate NICMOS 
PSF generated by Tinytim (Krist et al. 1998). Here we 
applied the same procedure for the NICMOS NIC1 F110W 
images (1.1 /im) and the WFPC2 optical images. In addi- 
tion we obtained nuclear WFPC1 or WFPC2 fluxes for a 
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number of galaxies from Peng & Ho (2001), see Tables 1 
and 2. These authors fitted a two dimensional model con- 
sisting of an analytical description for the bulge light plus 
an additional point source for the nucleus, convolved with 
an appropriate synthetic PSF. Their modeling for the un- 
resolved component is thus similar to ours. 

4.2. K' and L-band images 

We estimated the non-stellar fractions at K' and L band 
by fitting the PSF surface brightness profile to the nu- 
cleus of the galaxy, and assumed that the fitted PSF in- 
tegrated flux represents the unresolved component (see 
Alonso-Herrero et al. 1998 for more details). The FWHMs 
of the PSF in the K' and L band were measured from stan- 
dard star observations and from stars present in the images 
(only in K'). We also produced surface brightness profiles 
of the host galaxies to determine the amount of extended 
emission. 

In Table 4, column (2) gives the measured FWHM of the 
nuclear emission in the L band, column (6) is the fraction 
of unresolved emission in a 3"-diametcr aperture, and col- 
umn (7) is the nuclear morphology in the L band. As can 
be seen from this table the fraction of unresolved emission 
in the L band within a 3"-diameter aperture varies from 
16 to 100%, although in most cases is within 40 — 70%. 
McLcod & Rieke (1995) also obtained estimates of the un- 
resolved fluxes in the K band. We find that our estimates 
agree with theirs to within the photometric and fitting 
errors of both works. 

4.3. Mid-infrared wavelengths and contribution from star 
formation to ISO fluxes 

In AGNs past ~ 3 /im we expect the contribution from 
the underlying old stellar population to the nuclear fluxes 
to be much reduced compared to that at shorter wave- 
lengths. However if the nuclear stellar population is domi- 
nated by a starburst there may be contamination at wave- 
lengths > 3 /im. We note that recently Krabbe et al. 
(2001) and Lira et al. (2001) have obtained ground-based 
N imaging of small samples of Seyfert galaxies and found 
that the mid-infrared 10 /im emission is dominated by a 
central unresolved source; however, a low surface bright- 
ness extended component powered by star formation might 
not be apparent in such data. 

Because of the relatively large ISO/ISOCAM aperture 
(diameter of 18", see Clavel et al. 2000) we need to deter- 
mine whether the 6.75 and 9.63 /im ISO fluxes are contam- 
inated by circumnuclear star formation. To this end we 
compare 6 cm radio measurements through a large (~ 15") 
aperture with those through a small (~ 0.3") aperture. 
The radio measurements are taken from Kukula et al. 
(1995) for the small aperture and from Rush, Malkan, & 
Edelson (1996) for the large aperture. The small aperture 
radio measurements of Seyfert galaxies are most likely to 
originate in the AGN whereas the large aperture fluxes will 
include contributions from both the AGN and extended 
star formation. Small 6 cm / / f ii ratios are therefore 
a warning that there may be a significant level of circum- 
nuclear star formation within the large aperture. 

In Table 5 we list Seyfert galaxies in the CfA sample 
with available ISO fluxes and fractions of nuclear to 15" 

2 Typically the apertures of ground-based data were of the order of 



radio emission of < 15%. For comparison in this table 
we also list, when available, the 10 /im ground-based 2 to 
ISO flux ratios, and whether there is evidence for extended 
emission in the L band and ISO images. The latter are 
from Clavel et al. (2000) and our own analysis of the im- 
ages. We have also included in this table NGC 5929 and 
NGC 5940 as they appear to have extended mid-infrared 
emission from the comparison of the ISO and ground- 
based 10 /im fluxes. 

Mkn 334, Mkn 841 and NGC 4051 have small L „////, K> , 

O.o ' 15 

ratios at 6 cm, but do not appear to be extended in ei- 
ther the ISO or the L-band images, and have mid-infrared 
ground to ISO ratios consistent with most of this emission 
stemming from a small region. NGC 4388 shows evidence 
for extended radio emission, but the ISO images reveal the 
presence of a bright point source on top of an extended dif- 
fuse component, so it is likely that most of the ISO fluxes 
are unresolved emission originating in the AGN. 

NGC 1144, NGC 3982, NGC 5033 and Mkn 266SW are 
cases where the ISO mid-infrared fluxes may contain a sig- 
nificant contribution from circumnuclear star formation. 
In NGC 1144 there are Hn regions a few arcseconds SW 
of the nucleus (see Figure 1 and also Hippelein 1989). In 
NGC 3982 there are H n regions within the 18" ISO aper- 
ture (Gonzalcz-Delgado et al. 1997), whereas in NGC 5033 
bright H n regions are present in the spiral arms to within 
5" of the nucleus (Evans et al. 1996). Mkn 266SW is part 
of an interacting system whose nuclei are separated by 
~ 10". Since we detect L-band emission from both nuclei 
(see Table 4) it is possible that the ISO fluxes also con- 
tain contributions from both nuclei. For these galaxies we 
list in Table 6 small aperture ground-based mid-infrared 
fluxes. If these are not available then the ISO fluxes are 
given as upper limits to the unresolved emission. Although 
there may still be some contribution from a compact star- 
burst even for the small aperture ground-based measure- 
ments, it will be energetically less important than that 
from extended star formation (Imanishi 2002). 

5. OBSERVED SPECTRAL ENERGY DISTRIBUTIONS 

In Tables 6 and 7, columns (2) through (9), we have 
summarized the optical and near- and mid-infrared un- 
resolved (non-stellar) fluxes for the CfA Seyfert 2s and 
Seyfert Is, respectively. Errors in all the unresolved fluxes 
are on average ±30%. We also indicate in these two tables, 
when appropriate, the upper limits to the nuclear fluxes. 
Because of the large aperture of the ISO/ISOPHOT 16 /im 
measurements, we show all the fluxes at this wavelength as 
upper limits as they may be contaminated to some degree 
by emission from the host galaxy. 

In Figure 2 we present the non-stellar SEDs of 
NGC 1068 and NGC 4151, where the 1.1 through 2.2 /mi 
data points are from 7LST/NICMOS, and the 3.8 and 
4.8 /im are ground-based data, all taken from Alonso- 
Herrero et al. (2001) and references therein. We also show 
the ISO spectrum of the AGN component of NGC 1068 (Le 
Floc'h et al. 2001), as well as the JSO/ISOPHOT 16/zm 
data points from Perez Garcia & Rodriguez Espinosa 
(2001). The additional mid-infrared wavelength points are 
ground-based observations through small (6 — 8.5") aper- 
tures (Rieke & Low 1972; 1975a; 1975b) and diffraction 

.-8.5". 
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limited measurements of the central emission (Alloin et 
al. 2000 for NGC 1068 and Radomski et al. 2002 for 
NGC 4151). These two galaxies provide a guide to the 
behavior of the other galaxies in the sample, for which 
less complete data are available. From 10 to 34/tm, the 
SEDs are very similar, and both are flat in vf v . From 1 
to 10/tm, the SED of NGC 1068 has the appearance of a 
very strongly reddened version of the SED of NGC 4151. 
Qualitatively, this difference is as predicted by the unified 
model. 

The non-stellar optical through mid-infrared (~ 16 /im) 
SEDs of the remaining galaxies are presented in Fig- 
ures 3, 4 and 5 for the CfA Seyfert 1 — 1.5s, Seyfert 
1.8 — 1.9s and Seyfert 2s, respectively. When plotting the 
SEDs and fitting the spectral indices (next section) the 
ground-based iV-band fluxes (derived to give monochro- 
matic fluxes of Rayleigh- Jeans spectra) need to be cor- 
rected to v f v = constant by multiplying them by a factor 
of 1.22 (see also Edelson et al. 1987). 

Variability in the near-infrared and optical region (e.g., 
Clavel, Wamsteker, & Glass 1989) is a concern when as- 
sembling SEDs of AGN from multiple data sources. Most 
of these data were observed within a time period of a 
couple years: the ISO lifetime was Nov. 1995 through 
May 1998, and the NICMOS lifetime prior the cryocooler 
era was Feb. 1997 through Dec. 1998. Most of the WFPC2 
images were observed within a year of these operational 
lifetimes. A number of the CfA galaxies have duplicate 
NICMOS observations. NGC 5273, NGC 5033, NGC 5347, 
UGC 12138, and UM 146 displayed variable levels (over a 
time scale of a few months) of the unresolved emission 
at 1.6/im but Mkn 573, NGC 5252 and Mkn 471 did not 
(Quillen et al. 2000). Typical levels of variations in the 
near-infrared (at 1.6 /im) fluxes are up to ~ 20% (Quillen 
et al. 2000). This level of variation should be within the 
uncertainties of the unresolved flux estimates in the near- 
infrared, so we do not expect that our SEDs arc strongly 
affected by intrinsic variability. 

A common parametrization to describe the non-stellar 
SEDs of AGN is to fit a spectral index (a) so that the 
observed fluxes can be expressed as power laws, /„ oc v~ a . 
For a significant fraction of the Seyfert 2s we only have 
upper limits to the unresolved emission - at short wave- 
lengths if a point source was not detected, and at longer 
wavelengths if there is evidence for stellar emission. We 
use the package asurv (La Valley, Isobe, & Feigelson 1992) 
to fit the SEDs because it allows us to include upper lim- 
its. We have fit spectral indices in the optical to infrared 
(up to 16 /im) spectral range (a op t-iR) and in the infrared 
(1 — 16 /im) spectral range alone (am). The spectral in- 
dices are presented in the last two columns of Tables 6 and 
7. In Figure 6 we show the distributions of the fitted spec- 
tral indices in terms of the Seyfert type as derived from 
optical spectroscopy. 

5.1. The effects of extended foreground extinction in the 

host galaxy 

Before we discuss the observed SEDs of the CfA Seyferts 
and compare them with torus models, we must consider 
the effects of dust in the host galaxy. Most of the Seyfert 
galaxies in the CfA sample have dust near their nuclei 
(Martini & Pogge 1999; Pogge & Martini 2002). Gener- 
ally the extinction is not strong enough to be prominent in 



the NICMOS images (see Martini & Pogge 1999), which 
sets a limit on the total extinction that can be present on 
scales resolvable by HST. As can be seen from their color 
maps, the dust is generally located in front of the nucleus 
(see also Regan & Mulchaey 1999). We can therefore con- 
sider it as a purely foreground medium. Martini & Pogge 
(1999) used the V — H colors of the Seyfert 2s in the CfA 
sample to derive moderate levels of optical extinction, up 
to Ay ~ 5 mag. 

To illustrate the effects of dust in the host galaxy, we 
took SEDs produced by the tapered torus+conc mod- 
els (see Section 6 and Table 8) from Efstathiou et al. 
(1995) and inferred their appearance as seen through a 
foreground dust screen. We consider three different view- 
ing angles to the AGN: 9 V = 0° (polar view of the AGN), 
6 V = 30° (intermediate view), and 8 V = 90° (equatorial, 
through the torus, view of the AGN). The intermediate 
viewing angle has been chosen to coincide with the half 
opening angle of this particular torus model, which corre- 
sponds to the transition between type 1 and type 2 ob- 
jects. For the foreground extinction we used values of 
Ay = 0.4, 1 and 5 mag, and the Rieke & Lebofsky (1985) 
extinction curve. The resulting spectra are shown in Fig- 
ure 7. The reddened spectrum from a nearly face-on (po- 
lar view) torus (underlying Seyfert 1, 9 V = 0°) will look 
similar to the intrinsic spectrum at a viewing angle closer 
to a true Seyfert 2 (9 V = 90°) given sufficient levels of 
foreground extinction (Ay > 5 mag). However, moder- 
ate amounts of extinction will not have a significant effect 
on the nuclear fluxes at wavelengths greater than approx- 
imately 2 /im. In addition, if a foreground screen of dust 
that is generally thick enough to modify the infrared SEDs 
has small holes, the view of the torus through these holes 
will dominate the flux in the near infrared and give an 
impression of a low level of extinction overall. 

The torus models place the source of the near-infrared 
emission close to the broad-line region, and of a diame- 
ter of order 1 pc (compared with the size of about 0.01 
pc for the BLR). This general geometry is supported by 
observations of near-infrared variability with timescales of 
order a year (e.g., Rieke & Lebofsky 1979; Lebofsky & 
Rieke 1980; Glass 1998 and references therein). To ob- 
scure the entire near-infrared-emitting region with a fore- 
ground cloud that nonetheless lets the broad emission lines 
through would take a highly contrived configuration for the 
obscuring cloud: a small hole aligned on the BLR with no 
other holes elsewhere where near-infrared light could es- 
cape in our direction. Thus, it is improbable that dust 
in the host galaxy could significantly modify the infrared 
SED and still allow broad lines to be detected in optical 
spectra. 

5.2. Seyfert 1 - 1.5s 

The most remarkable aspect of the optical to mid- 
infrared non-stellar SEDs for the type 1 — 1.5 galaxies (see 
Figure 3) is the high degree of uniformity, with a feature- 
less continuum that is nearly flat in vf v , similar to that of 
NGC 4151. Overall, the newly derived SEDs are not dra- 
matically different from those inferred in previous works, 
since the dominance of the AGN emission over the stel- 
lar emission in the nuclear regions makes detailed image 
deconvolution unimportant to derive true nuclear fluxes. 
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However if we compare our SEDs on a case-by-case ba- 
sis with those in Edclson et al. (1987), it is clear that 
their Seyfert 1 SEDs between 1 and 16 /im tend to ap- 
pear bluer in v /„ than ours. We believe the fluxes at 
shorter wavelengths in Edelson et al. (1987) may still con- 
tain some contribution from stellar emission (these authors 
used 8.5"-diameter apertures for their photometry). 

The average values for the optical-infrared and infrared 
spectral indices of the Seyfert 1 — 1.5s (22 galaxies) are 
Oopt-iR = 1-59 ± 0.30 and ckir = 1-48 ± 0.30, where we 
have excluded NGC 6104 and Mkn 993 (optically classi- 
fied as Seyfert 1.5 — 1.8 and Seyfert 1.5 — 2, respectively) 
as they clearly show SEDs more compatible with those 
of type 2 objects. For comparison Edelson et al. (1987) 
found a = 1.15 ± 0.29 between 2.2 /jm and 25 /im for the 
CfA sample, and Fadda et al. (1998) measured a me- 
dian spectral index between 1 /im and 20 ^m of a = 1.4 
for a heterogeneous sample of Seyfert Is. The fact that 
the infrared spectral indices are on average bluer than the 
optical-infrared indices is probably due to the presence of 
some amount of foreground dust in Seyfert Is. 

However, two type 1 galaxies, Mkn 590 and NGC 7603, 
have well-measured SEDs that depart significantly from 
the pattern established for the other type 1 objects. The 
SEDs of both start to drop from the near infrared towards 
the visible, relative to the fiat SEDs (in vf v ) that are typ- 
ical. We quantify this behavior from the nuclear fluxes 
measured at K' and L since they use identical observa- 
tional techniques and were obtained simultaneously (hence 
variability is not an issue) . Relative to the median K' — L 
for the type 1 nuclei, an extinction of order Ay = 10 mag 
is indicated for both galaxies. A similar extinction level 
would apply to the BLR, which would make the identi- 
fication of these galaxies as type 1 very improbable. We 
conclude that it is implausible that the differences in the 
SEDs for these two galaxies compared with the other type 

1 objects can be explained in terms of identical intrinsic 
SEDs modified by foreground extinction. 

5.3. Seyferts 2s 

In clear contrast with the CfA Seyfert 1 — 1.5s, the 
Seyfert 2s show a variety of spectral shapes (see Figure 5, 
and also Figure 6). All the Seyfert 2 galaxies have SEDs 
that rise steeply at increasing wavelengths (chir, > 2.6), 
qualitatively similar to the SED of NGC 1068. This be- 
havior is roughly consistent with the unified model predic- 
tions for a Seyfert 2 nucleus, that is, a heavily obscured 
version of the typical Seyfert 1 SED. 

The Seyfert 2 nuclei in the CfA sample have near- 
infrared colors and near-to-mid infrared colors in general 
redder than seen in the compilation by Fadda et al. (1998), 
but similar to the Seyfert 2s observed by Alonso-Herrero et 
al. (2001). One possibility is that the Seyfert 2 photome- 
try compiled by Fadda et al. (1998) was still contaminated 
with light from the galaxy for the short wavelengths, as the 
determination of the unresolved component is usually lim- 
ited by the seeing conditions when using ground-based ob- 
servations (see Alonso-Herrero et al. 1996). The infrared 
spectral index of the median SED for a type 2 object in 
Fadda et al. (1998) is approximately aiR = 2.2 between 1 
and 20 /im. However from Figures 5 and 6 it is clear that 
it is not appropriate to define a 'median' SED for type 

2 objects, as it is precisely the variety of spectral shapes 



that needs to be reproduced by the torus models, as we 
will discuss in Section 7. 

Some of the type 2 nuclei have substantial amounts of 
extended foreground dust associated with their host galax- 
ies. For instance NGC 4388 and NGC 5929 have measured 
NLR extinctions of Ay ~ 4 mag (see Veilleux et al. 1997 
and Rhee & Larkin 2002, respectively). The effect of cor- 
recting the observed SEDs for the indicated amounts of 
foreground extinction is small. Foreground dust at the 
level seen in these galaxies can strongly affect the optical 
classification of the nuclear activity if it preferentially ob- 
scures the broad-line region (as proposed by Maiolino & 
Rickc 1995). However, this amount of dust does not seem 
capable of influencing the infrared SEDs sufficiently to ac- 
count for the systematic differences seen between Seyfert 
1 and Seyfert 2 SEDs. 

5.4. Seyfert 1.8 - 1.9s 

UGC 12138, NGC 3786, NGC 5033, and NGC 5273 
(Figure 4) are optically classified as Seyfert 1.8 — 1.9s, but 
have nearly flat SEDs, similar to those of Seyfert 1 — 1.5s 
(see Figure 3). In fact their infrared spectral indices fall 
well within the values found for the CfA Seyfert 1 — 1.5s. 
The extinction maps of Martini & Pogge (1999), based on 
V — H colors, show that all four of these galaxies are free of 
heavy, extended obscuring foreground clouds in their host 
galaxies. 

We exclude from classification as type 1.9 those galax- 
ies where the broad line component is only seen in the in- 
frared, e.g., Pa/3. The intermediate-class galaxies UM 146, 
Mkn 334, Mkn 471, NGC 4258, NGC 5252, and NGC 5674 
(Figure 4), display intermediate SEDs (infrared spectral 
indices in the range ckir, ~ 1.8 — 2.6) that do not drop to- 
ward short wavelengths as steeply as NGC 1068. In most 
of these galaxies, unresolved nuclear emission is measured 
from optical and/or near-infrared HST images. In all of 
these cases, the V — H images of Martini & Pogge (1999) 
suggest a possibility of extended obscuration from the host 
galaxy, although the bright and blue nucleus of Mkn 471 
indicates that it is probably clear of the nearby clouds. 

One possibility for the Syl.8 — 1.9 intermediate SEDs is 
that they are those of a true Seyfert 1 galaxy seen through 
some foreground dust. However, an explanation of the 
steep infrared SEDs in terms of foreground dust is un- 
likely in the case of NGC 4258. We have used the same 
approach described under type 1 objects, estimating the 
extinction needed to redden a median type 1 infrared SED 
to agree with the K' — L colors, and find for this galaxy 
a level equivalent to Ay > 10, making observation of any 
broad-line component very problematic. A similar conclu- 
sion holds for NGC 7479; its K' - L is so red that the 
nucleus is not detected at the shorter wavelengths, but 
the NICMOS upper limits demonstrate that it is heavily 
obscured (if the intrinsic SED is Seyfert- 1-like). 

We therefore find four cases among the type 1 and 
intermediate-class Seyfert galaxies where broad line com- 
ponents are seen despite an infrared SED that is either 
intrinsically red or is heavily reddened. The conclusion 
that broad line components can escape from systems with 
such SEDs is reinforced by the detection by Ho et al. 
(1997) of broad wings on Ha in NGC 1068, NGC 3982, 
and NGC 4388. For many of this total of seven objects, 
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under the unified model the "extinction" we see in the in- 
frared SED is likely to be a property of the torus itself, 
rather than the result of a foreground obscuring cloud. 

The Scyferts with intermediate SEDs are the most diffi- 
cult for the disk and torus models to explain. The models 
are extremely sensitive to the viewing angle. Only for a 
very small range of orientation angles, just grazing the 
surface of the disk or torus, are intermediate SEDs pre- 
dicted. The models would also predict that the broad 
emission lines would be obscured at virtually all viewing 
angles where the near-infrared emission region is obscured. 
Yet we find that 10 — 20% of the galaxies with detectable 
broad lines in the optical have red infrared SEDs. That is, 
a strict dichotomy of SEDs correlated with emission-line 
properties is predicted, but is not in accordance with our 
measurements. 

5.5. Spectral indices and the inclination of the host 

galaxy 

There appears to be a deficit of Seyfert 1 — 1.5 galax- 
ies in edge-on galaxies (e.g., Keel 1980; Lawrence & Elvis 
1982; Maiolino & Rieke 1995; Schmitt et al. 2001) in most 
of the optically selected samples. In particular, Maiolino 
& Rieke (1995) concluded that the CfA sample has a de- 
ficiency of Seyfert Is in edge-on galaxies, caused by an 
extended torus (scale of a 100 pc) coplanar with the host 
galaxy that obscures the BLR. They predict that the spec- 
trum of a Seyfert 1 galaxy will look Seyfert 2-like in many 
of these cases. 

The behavior of the infrared SEDs provides an indepen- 
dent way to estimate whether AGNs are of type 1 or 2. 
Where the nuclei may be obscured by dust in the host 
galaxy, the SEDs are much less affected than the opti- 
cal emission lines. Using the SED shapes, we can check 
whether the portion of intrinsically type 1 galaxies is sim- 
ilar in edge-on systems to the portion found in optically 
selected samples in face-on galaxies. 

The axial ratios (b/a) for the extended CfA sample of 
Seyfert galaxies are taken from McLeod & Rieke (1995) 
and Ho et al. (1997). In Figure 8 we show the spectral 
indices as a function of the host galaxy inclination, b/a. 
The galaxies have been divided in Seyfert Is (1 and 1.5 
types) and Seyfert 2s (1.8, 1.9 and 2 types) according to 
their optical classification. We also mark in this figure for 
reference an axial ratio of b/a — 0.6 (or inclination of the 
host galaxy of i — 53°). Maiolino & Rieke (1995) found 
that pure Seyfert 1 galaxies with axial ratios of b/a < 0.6 
will be most likely classified as Seyfert 1.8 or 1.9. This 
result is confirmed by Figure 8. 

If we use the infrared SED to divide galaxies into prob- 
able type 1 and type 2, we can see that in terms of the 
spectral index «ir there are approximately equal numbers 
of type 1 (aiR = 1.52 ± 0.38, 6 galaxies) and type 2 (5 
galaxies) in the highly inclined galaxies in the extended 
CfA sample (Figure 8, right panel). These numbers are 
consistent those for the face-on members of the sample. 
This result confirms the suggestion by Maiolino & Rieke 
(1995) that Seyfert 1 nuclei hosted in highly inclined galax- 
ies may be misclassified as type 2 objects. For example, 
studies based on optical spectroscopy and that ignore in- 
clination effects will overestimate the portion of type 2 
nuclei. 



6. MODEL PREDICTIONS 

For the following discussion, we assume the validity of 
a unified model, that is we assume that all Seyfert galax- 
ies have similar dusty tori in their nuclei. We will try to 
find a geometry for these tori that is consistent with our 
observations. 

The broadness of the spectrum in the 1 — 20 /im region 
and the nature of the 9.7 fim silicate feature (in absorption 
and emission) have been a challenge for the torus mod- 
els. Compact torus models (Pier & Krolik 1993) often 
predict spectra that are not sufficiently broad (they lack 
emission in the near-infrared) or they emit too strongly 
in the 9.7 fxm silicate feature for polar orientations. This 
issue has resulted in a suite of more extended disk models: 
the flared disk (Granato et al. 1997), the tapered disk 
(Efstathiou & Rowan- Robinson 1995), and the warped 
disk (Sanders et al. 1989). Excess emission in the near- 
infrared in NGC 1068 prompted Pier & Krolik (1993) to 
suggest that there was an additional clumpy dust compo- 
nent above the torus. This two component torus+cone 
model was also used by Efstathiou et al. (1995) to match 
the broad spectrum of NGC 1068. 

To suppress the silicate emission in type 1 objects, pre- 
vious models forced the torus geometry to have extremely 
steep inner regions so that the region directly illuminated 
by UV radiation could only be viewed at high inclination 
angles. Since the models were also driven by the belief 
that there was strong silicate absorption in the Seyfert 2s, 
the outer disks were arranged so that at high inclination 
angles the inner hotter regions were only viewed through 
colder material (e.g., Pier & Krolik 1993; Efstathiou & 
Rowan-Robinson 1995; Granato et al. 1997). Also Pier & 
Krolik (1993) models assumed high optical depths to be 
consistent with the high hard X-ray column densities often 
measured in Seyfert 2 galaxies. 

Clavel et al. (2000) have used ISO spectroscopy 
(through a 24" x 24" aperture) to show that Seyfert 1 
galaxies have little silicate emission, and a broad smooth 
infrared continuum. However, the Seyfert 2s displayed 
PAH emission and lower levels of continuum. In fact re- 
cent ISO and ground-based mid-infrared spectroscopy of 
a small number of nearby AGN show that the silicate ab- 
sorption feature in some Seyfert 2s might be shallower than 
previously inferred (see e.g., Laurent et al. 2000, Imanishi 
& Ueno 2000 and Le Floc'h et al. 2001) because measure- 
ments of the nearby continuum were contaminated with 
PAH features originating in the host galaxy. The mid- 
infrared emission is probably slightly extended in many 
Seyfert 2s (Maiolino et al. 1995) across the region covered 
by the ISO spectra, so it is likely that there is some PAH 
contribution from the surrounding galaxy in a number of 
cases. The nuclear spectrum itself might contain little or 
no PAH emission. For example, little silicate absorption is 
seen in the ISO/SWS nuclear spectrum of NGC 1068 and 
the continuum is nearly featureless (Le Floc'h al. 2000, 
see also our Figure 2, left panel). 

In addition to the behavior of silicate feature, torus mod- 
els can be constrained by the shapes of the SEDs for a 
complete sample of galaxies. To do so, we use the tapered 
disk model of Efstathiou & Rowan- Robinson (1995). In 
this geometry, the height of the disk increases with the 
radial distance, but tapers off to a constant height in the 
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outer part. The main parameters of the model arc the half 
opening angle of the torus O c , the equatorial UV optical 
depth t\jv, the dust sublimation temperature T, the ratio 
between inner and outer radii ri/r 2 , and the ratio between 
the scale height and the outer radius hjr 2 . Table 8 summa- 
rizes the parameters of the models used here. In addition 
to the torus models, we will also consider the composite 
torus+cone model used by Efstathiou et al. (1995) to fit 
the SED of NGC 1068. The viewing angles are defined as 
9 V = 0° for the polar view of the AGN, and V = 90° for the 
equatorial view of the AGN (note that in the Efstathiou 
& Rowan-Robinson 1995 models the authors measured the 
angles from the equator). 

Because of their high equatorial opacities, the torus 
models presented here (but also those of Pier & Krolik 
1993; Granato & Dancse 1994; Granato et al. 1997) gen- 
erally predict the following: (i) for viewing angles 9 V < 9 C 
an SED more or less flat in vj v depending on the model 
input parameters, and a direct view of the nuclear optical 
emission; (ii) for viewing angles 9 V > O c , a steep near and 
mid-infrared SED becoming steeper as the viewing angle 
increases, with no direct view of the nucleus; (iii) an inter- 
mediate SED as well as some nuclear optical emission for 
only a very small range of possible viewing angles around 
8 V ~ C . A dichotomy of spectral shapes is thus predicted, 
and few, if any, Seyferts are expected to have an interme- 
diate SED. This behavior is shown graphically in Figure 9 
where the predicted infrared spectral index («ir) is plot- 
ted against the cosine of the viewing angle for some of the 
models described in Table 8. 

7. COMPARISON WITH OBSERVATIONS 

The observed distribution of q;ir (see Figure 6) in the 
extended CfA sample clearly lacks the dichotomy of spec- 
tral indices predicted by the models. As discussed in Sec- 
tion 5, all the Seyfert 1 — 1.5s and a large fraction of the 
Seyfcrt 1.8 — 1.9s show infrared spectral indices clustered 
near am, — 1.5 ± 0.3. Type 2 objects on the other hand 
show a variety of steeper spectral indices, unlike the model 
predictions. Moreover the large opacity torus models with- 
out the extended cone component predict steeper SEDs 
(«ir > 4) than observed in the CfA Seyferts for suffi- 
ciently large viewing angles. Below, we discuss a number 
of possibilities to reconcile the observed SEDs with the 
model predictions. 

7.1. Torus models with low opacities 

Figure 9 shows that even the model with an equatorial 
optical depth of tuv = 250 mag (Ay = 50 mag) produces 
SEDs with aiR > 4 for viewing angles 9 V > 60° . The other 
torus models considered here with tuv > 250 mag predict 
even steeper SEDs near an equatorial view of the AGN. 
Granato et al. (1997) and Fadda et al. (1998) used simi- 
lar arguments based on the observed near and mid-infrared 
colors of Seyfert galaxies to favor low opacity (equatorial 
optical depths Ay — 10 — 30 mag) for extended (a few 
hundred parsec) torus models. 

Traditionally low opacity torus models have been ruled 
out because they tend to produce the 9.7 /im silicate 
feature in emission for viewing angles of less than the 
half opening angle of the torus (that is, type 1 objects). 
Ground-based mid-infrared spectroscopy of Seyfert Is did 
not show evidence for this feature in emission (Roche et 



al. 1991). Recently, however, Clavel et al. (2000) have 
argued that the average ISO mid-infrared spectrum of a 
sample of Seyfert 1 galaxies displays some evidence for the 
presence of the silicate feature in emission. Spoon et al. 
(2002) on the other hand do not find evidence for this fea- 
ture. Additional observations are needed to resolve this 
issue. 

7.2. Torus models with an extended (optically thin) 
component 

The extended conical component proposed in Efstathiou 
et al. (1995) was assumed to be optically thin and thus 
increases the continuum level around 3 — 5 /im, causing 
the infrared spectral index to be shallower even for high 
viewing angles (see Figure 9). Indeed, NGC 1068 shows 
emission by dust from the wall of the cavity, i.e., perpen- 
dicular to the axis of the torus and in an extended com- 
ponent (100 pc) coincident with the radio emission and 
ionization cone (see Bock et al. 1998). The Torus+Cone 
model of Efstathiou et al. (1995) seems to reproduce well 
the observed distribution of ojir in the type 2 objects in 
the CfA sample, although it will be necessary to increase 
the opening angle of the torus to match the observed rel- 
ative numbers of type 1 and type 2 objects in the CfA 
sample. This model also reproduces well the observed in- 
frared spectral indices of Seyfert 1 — 1.5s. 

This model still fails to solve the problem of the inter- 
mediate Seyferts, as infrared spectral indices in the range 
c*ir = 1.5 ± 0.1 — 2.5 ± 0.1 are extremely unlikely to 
occur. This discrepancy can be solved with fine tuning 
of the models. For instance, a torus+cone model with 
a lower equatorial optical depth would make the transi- 
tion of infrared spectral indices between type 1 and type 2 
smoother, thus narrowing the range of 'unlikely' spectral 
indices. Another way of producing intermediate SEDs is 
having a range in the visibility of the scattering region. 
This could be done with a range in the h/r 2 parameter 
or if the scattering region is physically located at differ- 
ent distances from the nucleus. A more visible scattering 
region would produce a stronger 1 — 2 /im continuum and 
therefore shallower SEDs. This could account for some of 
the intermediate Seyferts (Figure 4). 

7.3. Clumpy torus models 

Nenkova, Ivezic, & Elitzur (2002) have recently put for- 
ward a more realistic model for the emission by dust in 
AGNs following the initial suggestion by Pier & Krolik 
(1993) that the dust in the torus must be in clumps to pro- 
tect the grains. The parameters of this clumpy torus are, 
apart from the usual geometry parameters of the torus, 
the radial distribution of clouds, the number of clouds and 
the optical depth of the individual clouds. In their model 
a cloud is heated by radiation from both the AGN and 
other surrounding clouds. The main result from this pa- 
per is that this model can reproduce simultaneously the 
broadness of the infrared SEDs and suppress the silicate 
feature emission with a few clouds (5 — 10) with optical 
depths Ty > 60. 

Nenkova et al. (2002) only showed results from the lim- 
iting viewing angles of 8y = 0° and 9y = 90°, but from 
their figure 3 we can see that the equatorial view of the 
AGN in that particular model produces an SED consistent 
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with the steepest SEDs in the CfA sample. It remains to 
be seen if their model could reproduce the variety of ob- 
served SEDs measured for type 2 objects in the CfA sam- 
ple. One concern in their comparison with observations, 
although they will provide a more detailed modelling in 
the near future, is the small number of type-2 SEDs fitted, 
and the possibility that some of the examples are contam- 
inated by stellar processes. For example, for Circinus, one 
of the Seyfert 2 galaxies used in their modelling, Ruiz et 
al. (2001) demonstrate that a dominant starburst compo- 
nent is present at wavelengths longer than approximately 
25 /urn. 

7.4. Dust composition 

The lack of deep silicate absorption features in Seyfert 
2s and other arguments based on anomalous Ay /Ah ratios 
in AGNs have prompted the idea that the composition and 
properties of dust in AGNs may be different from those of 
Galactic dust (see e.g., Laor & Draine 1993; Maiolino et 
al. 2001; Maiolino, Marconi, & Oliva 2001). 

Specifically, the dust around AGNs is proposed to have 
a large portion of large grains (Maiolino et al. 2001: 
Maiolino, Marconi, & Oliva 2001). In that case, the ex- 
tinction curve would be flatter than the standard Galactic 
extinction curve. A model with a higher fraction of large 
grains (or fewer silicates) means we can have a flat type 1 
spectrum (with weak or no silicate emission feature) with a 
lower optical depth, as well as shallower SEDs for viewing 
angles 8 V > O c (that is, type 2s). Fadda et al. (1998) have 
shown the effects on near infrared and mid infrared colors 
using models with fewer silicate grains than the standard 
Galactic dust composition. Such models would be in bet- 
ter agreement with our findings for the extended CfA sam- 
ple. 

8. SUMMARY 

We define an extended CfA sample of Seyfert galaxies 
composed of the original CfA Seyferts and those galax- 
ies initially classified as LINERs, but later reclassified as 
Seyferts from higher S/N optical spectra. For this sample, 
we present high resolution IRTF K 1 (2.2 /im) and L-band 
(3.5 /im) images and compile HST, ISO and ground-based 
observations in the spectral region 0.4—16 /im. These data 
have been used to determine the unresolved nuclear emis- 
sion with improved accuracy over any previous studies. 

The resulting spectral energy distributions are used to 
test and constrain unified models for the dusty tori in 
the active nuclei. We find that the fitted infrared spec- 
tral indices (/„ oc v~ am ) between 1 and 16 /mi are corre- 
lated with the Seyfert type in the sense that steeper SEDs 
tend to be found in Seyfert 2s and flatter ones in Seyfert 
1 — 1.5s. The majority of galaxies optically classified as 
Seyfert types 1.8 and 1.9 display infrared spectral indices 
and SEDs either similar to those of type 1 objects or inter- 
mediate between Seyfert Is and Seyfert 2s. The interme- 
diate SEDs of many Seyfert 1.8 — 1.9s may be consistent 



with the presence of a pure Seyfert 1 viewed through a 
moderate amount (Ay <^ 5 mag) of foreground galaxy ex- 
tinction. 

However, the abrupt division in SED shape correlated 
with emission-line properties that would be predicted by 
an optically thick torus model — flat for type Is and steep 
for type 2s — is not observed. Discrepancies with this 
model include the wide range of spectral indices observed 
in type 2 objects, the lack of extremely steep SEDs, and 
the large number of objects with intermediate properties. 
In addition, there are a surprisingly large number of galax- 
ies with broad-line components and steep SEDs, a com- 
bination that conventional torus models would predict is 
very rare. 

These results, combined with the lack of deep 9.7 /mi sil- 
icate absorption feature in Seyfert 2s and the possible pres- 
ence of silicate emission in Seyfert Is (Clavel et al. 2000), 
suggest that the region of infrared emission might not be as 
optically thick as previously considered. Although X-ray 
and maser observations show that high opacity material 
does exist near the AGN, it is not clear whether this ma- 
terial is responsible for the infrared emission. In a recent 
paper, Risaliti, Elvis, & Nicastro (2002) have shown that 
the observed fast variability (a few months) of hard X-ray 
column densities in Seyfert 2s is not consistent with the 
absorbing material being located in the standard parsec 
scale torus. Since this material may be closer to the black 
hole than the region responsible for the infrared emission 
in Seyfert galaxies, the hard X-ray column densities might 
not be useful for constraining the torus optical depth. 

Other possible explanations for the observed SEDs in- 
clude an optically thick torus with an additional optically 
thin infrared component (such as from a cone of dust 
clouds: Efstathiou et al. 1995) to account for the lack 
of very steep SEDs and the lack of deep silicate absorp- 
tion features. Another promising approach is the clumpy 
torus model of Nenkova et al. (2002). 
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FIGURE CAPTIONS 

Figure 1. — Contours plots of selected galaxies with extended emission in the L-band (right panels). We show the 
-ftT'-band contour plots for comparison (left panels). The orientation of the images is south up, east to the left. 

Figure 2. — Left panel: Non-stellar SED of NGC 1068 in units of erg cm" 2 s -1 . The ISO spectrum of the AGN 
(nuclear) component of NGC 1068 is from Le Floc'h et al. (2001). The additional mid-infrared data points at 12.6, 17, 
22.5, 24.5, and 33.5 /im are ground-based small aperture observations from Rieke & Low (1975a); data points at 11.2 and 
20.5 /im are from Alloin et al. (2000). Right panel: Non-stellar SED of NGC 4151 in units of erg cm -2 s _1 . The additional 
mid-infrared data points at 10 and 18 fim (unresolved emission) are from Radomski et al. (2002), at 21 /im is from Rickc 
& Low (1972) and at 33.5 ^m from Rickc & Low (1975b). 

Figure 3. — Observed non-stellar SEDs (v fv) of optically classified type 1 (1 — 1.5) CfA Seyferts in units of ergcm~ 2 s _1 . 

Figure 4. — Observed non-stellar SEDs (v /„) of optically classified intermediate type (1.8 — 1.9) CfA Seyferts in units 
of erg cm~ 2 s _1 . 

Figure 5. — Observed non-stellar SEDs (v f v ) of optically classified type 2 CfA Seyferts in units of ergcm~ 2 s _1 . 

Figure 6. — Distributions of the fitted spectral indices (/„ oc v~ a ) in the optical-infrared wavelengths (a pt-iR for the 
spectral range ~ 0.4 — 16 fim, left panels) and infrared wavelengths (am for the spectral range 1 — 16 ^m, right panels). 
The very top panels are the 51 Seyferts in the 'extended' CfA sample with measured spectral indices. The panels below 
are the same sample broken up in different Seyfert types as derived from optical spectroscopy. 

Figure 7. — Effects of foreground extinction (Ay = 0.4, land 5 mag) on the observed SED of an AGN, using Rickc 
& Lebofsky (1985) extinction law. The AGN SEDs are outputs of the Torus+Cone model from Efstathiou & Rowan- 
Robinson 1995 (see Section 6 and Table 8 for details) for three different viewing angles. 8 V = 0° is the polar view of the 
AGN, 6 V = 30° is an intermediate view (corresponding to the half opening angle of the torus, O c ), and 6 V — 90° is the 
equatorial view (through the torus). 

Figure 8. — Fitted optical-infrared (left panel) and infrared (right panel) spectral indices (see Section 5) as a function 
of the inclination of the host galaxy, where b/a = is an edge-on galaxy, and b/a — 1 is a face-on galaxy. The arrows 
at the top of the plots indicate the values of b/a for the remaining 7 galaxies with no measured spectral indices — dotted 
lines for the type 2s (green in the on-line version) and solid lines for the type Is (red in the on-line version). The dashed 
line indicates a host galaxy axial ratio of b/a = 0.6 or an inclination of i > 53° (see discussion in text). The galaxies are 
separated between Seyfert Is (1 — 1.5) shown as open circles (red in the on-line version) and Seyfert 2s (1.8 — 2) shown 
as filled star symbols (green in the on-line version), using the classification from optical spectroscopy. 

Figure 9. — Predicted infrared spectral indices (f v oc v~ aiIi ) in the 1 — 16 /im spectral range as a function of the 
cosine of the viewing angle to the torus (6 V = 0° polar view of the AGN, 8 V = 90° equatorial view of the AGN) using 
Efstathiou & Rowan-Robinson (1995) and Efstathiou et al. (1995) torus models. Left panel: results for models with a 
torus half opening angle 9 C = 45°, and two different equatorial UV optical depths: tuv = 500 mag (filled triangles), 
t\jv = 250 mag (filled stars), Models Torus 2 and Torus 3 respectively (see Table 8). Right panel: results for a torus model 
with a half opening angle O c = 30°, and an equatorial UV optical depth of tuv = 1200 mag with the cone component 
(Model Torus+Cone in Table 8, open circles) and without the cone component (Model Torus 4 in Table 8, filled circles). 
In both diagrams the shaded areas represent the observed ranges of fitted am for Seyfert 1 — 1.5s (darkest shaded area), 
Seyferts 1.8 — 1.9s (intermediate shaded area), and Seyfert 2s (lightest shaded area) in the CfA sample. 
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Table 1 

Summary of high resolution optical, near-infrared and mid-infrared observations for the CfA 

Seyfert galaxies. 



Galaxy 


Typo 


WFPC1/2 


Rcf 


NICMOS 


Rcf 


KL 


Mid-IR 


Rcf 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Mkn 334 


Syl.8 




1 


F110W, F160W 


2, 4 


2 


6.75, 9.63, 10, 16Mm 


5, 9, 6 


Mkn 335 


Syl.O 


F785LP 


1 






2 




5, 6 


UGC 524/0048+29 


Syl.O 












10/J.m 


9 


I Zw 1 


Syl.O 


F555W 


1 






10 


5.9, 7.7, 10, 16/J.m 


7, 10, 6 


Mkn 993 


Syl.5 - 2.0 


F606W 


1 






2 


10, 16jU.m 


11, 6 


Mkn 573 


Sy2.0 (non-HBLR) 


F708W 


2 


F110W, F160W 


2, 4 


3 


10. 16jt*m 


11, 6 


UM 146/UGC 1395 


Syl.9 


F606W 


1 


F110W, F160W 


2, 4 


2 


6.75, 10, 16/J.m 


2, 9. 6 


Mkn 590/NGC 863 


Syl.O 


F785LP 


1 






2 


6.75, 9.63, 16/J.m 


5, 6 


NGC 1068 


Sy2.0/Syl.8 (HBLR) 






F110W, F160W 


3 


3 


5 — 16, 16/^m 


13,6 










F222M 










NGC 1144 


Sy2.0 (non-HBLR) 


F606W 


2 


F110W, F160W 


2, 4 


2 


6.75, 10, le^m 


2, 11, 6 


Mkn 1243/NGC 3080 


Syl.O 












10. 16jum 


9, 6 


NGC 3227 


Syl.5 


F547M 


1 


F160W, F222M 


3 


2 


6.75. 9.63, 16Mm 


5, 6 


NGC 3362 


Sy2.0 (non-HBLR) 


F606W 


2 


F110W, F160W 


2, 4 


2 


10. 16 M m 


9, 6 


UGC 6100 


Sy2.0 (non-HBLR) 


F606W 


2 


F110W, F160W 


2, 4 


2 


6.75, 9.63, 10, 16^m 


5, 9, 6 


NGC 3516 


Syl.5 


F547M 


1 


F160W 


4 


14 


6.75, 9.63, 16/j.m 


5, 6 


NGC 3786/Mkn 744 


Syl.8 




1 


F110W, F160W 


2, 4 


2 


10. 16jum 


9, 6 


NGC 3982 


Sy2.0/Syl.9 (non-HBLR) 


F606W 


2 


F110W, F160W 


2, 4 


2 


6.75, 9.63, 10, 16^111 


5, 9. 6 


NGC 4051 


Syl.O 


F547M 


1 






2 


6.75, 9.63, 10, 16/im 


5, 10, 6 


NGC 4151 


Syl.5 


F547M 


1 


F110W, F160W 


3 


3 


5.9, 7.7, 10, 16/xm 


7, 10, 6 








F222M 








NGC 4235 


Syl.O 


F547M 


1 


F160W 


4 


2 


10, 16(im 


9, 6 


NGC 4253/Mkn 766 


Syl.5 


F785LP 


1 


F160W 


4 


2 


6.75, 9.63, 10, 16 M m 


5, 9. 6 


Mkn 205 


Syl.O 


F814W 


1 






10 


6.75. 10, lGfim. 


2, 16, 6 


NGC 4388 


Sy2.0/Syl.9 (HBLR) 






F110W, F160W 


2, 4 


8 


6.75, 9.63, 10, 16//m 


5, 17, 6 


NGC 4395 


Syl.0/Syl.8 


F450W 


1 


F160W 


4 


2 


10, 12/^m 


9, 15 


Mkn 231 


Syl.O 


F814W 


1 


F110W, F160W 


2, 4 


10 


5.9, 7.7, 10, 16/Ltm 


7, 10, 6 


NGC 5033 


Syl.9/Syl.5 


F547M 


2 


F110W, F160W 


2, 4 


2 


6.75, 9.63. 10. 16 M m 


5, 9. 6 


UGC 8621/1335 + 39 


Syl.8 












6.75. 9.63, le^m 


5, 6 


NGC 5252 


Syl.9 (HBLR) 


F547M 


2 


F110W, F160W 


2, 4 


3 


6.75, 9.63, 10, 16^m 


3, 9 


NGC 5256/Mkn 266SW 


Sy2.0 (non-HBLR) 






F110W, F160W 


2, 4 


2 


6.75, 9.63, 10, 16^m 


5, 9. 6 


NGC 5283/Mkn 270 


Sy2.0 (non-HBLR) 


F708W 


2 


F110W, F160W 


2, 4 


2 


6.75, 10, 16Mm 


2, 10, 6 


NGC 5273 


Syl.9/Syl.5 


F439W 


2 


F110W, F160W 


2, 4 


2 


6.75, 9.63, 10. 16jU.ni 


5, 9. 6 


Mkn 461 


Sy2.0 






F110W, F160W 


2, 4 


2 


6.75, 10/j-m 


2, 9 


NGC 5347 


Sy2.0 (non-HBLR) 






F110W, F160W 


2, 4 


2 


10/J.m 


9 


Mkn 279 


Syl.O 




1 






10 


6.75. 9.63, 16Mm 


5, 6 


NGC 5548 


Syl.5 


F785LP 


1 


F160W, F222M 


3 


2 


6.75, 9.63, 16/xm 


5, 6 


Mkn 471 


Syl.8 


F606W 


1 


F110W, F160W 


2, 4 


10 


10 M m 


10 


NGC 5674 


Syl.9 


F606W 


1 


F110W, F160W 


2, 4 


2 


6.75, 9.63, 10/xm 


5, 9 


Mkn 817 


Syl.5 


F785LP 


1 






12 


6.75, 9.63, le^m 


5, 6 


NGC 5695/Mkn 686 


Sy2.0 (non-HBLR) 


F606W 


2 


F110W, F160W 


2, 4 




10. 12fim 


9, 2 


Mkn 841 


Syl.5 


F785LP 


1 






2 


6.75, 9.63, 10. 16(im 


5, 9. 6 
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Table 1 — Continued 



Galaxy 


Typo 


WFPC1/2 


Ref 


NICMOS 


Ref 


K L 


Mid-IR 


Ref 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


O) 



NGC 5929 
NGC 5940 
Mkn 1513 
NGC 6104 
UGC 12138 

NGC 7469 

NGC 7603/Mkn 530 
NGC 7674/Mkn 533 
NGC 7682 



Sy2.0 (non-HBLR) 
Syl.O 
Syl.O 
Syl.5 - 1.8 
Syl.8 

Syl.O 

Syl.5/Syl.8 
Sy2.0 (HBLR) 
Sy2.0 (HBLR) 



F606W 


2 


F110W, F160W 


2, 4 


2 


6.75, 9.63, 10, 16^m 


5, 9, 6 


F606W 


1 


F160W 


4 


2 


6.75. 9.63, lOjum 


2, 9 










2 


6.75. 10, 14.3 M m 


2, 10 


F606W 


1 


F160W 


4 




10. 16 M m 


9, 6 




1 


F110W, F160W 


2, 4 


2 


6.75, 9.63, 10, 16 M m 


2, 9, 6 


F7S5LP 


1 


F110W, F160W 


3 


2 


6.75, 10, 16Aim 


8, 10, 6 






F222M 










F785LP 


1 






2 


6.75. 9.63. 10 M m 


5, 10 






F110W, F160W 


2, 4 


3 


6.75, 9.63, 10, 16/iAm 


5, 9, 6 


F606W 


2 


F110W, F160W 


2, 4 


2 


6.75, 10, 16^m 


2, 9, 6 



NOTE.— Column (2): Seyfert types taken from Osterbrock & Martel (1993) and Ho et al. (1997). NGC 6104 has been classified as a Seyfcrt 1.5 - 1.8 galaxy by Jimenez- 
Benito et al. (2000). Goodrich (1995) reports that NGC 7603/Mkn 530 and Mkn 993 (see also Tran, Osterbrock, Sz. Martel 1992) have undergone transitions between types 
Seyfert 1 and Seyfcrt 1.8 — 1.9. Wc also indicate in this column for the Seyfcrt 2s whether hidden broad line regions have been detected from spectropolarimctry (HBLR) 
or not (non-HBLR). Column (3): HST optical filter, unless the * symbol is used in which case the optical unresolved flux was estimated from ground-based observation (see 
Ho &6 Peng (2001) for more details). Column (4): references for the unresolved optical fluxes. Column (5): HST/NICM.OS filters. Column (6): references for the NICMOS 
fluxes. Column (7): references for the K L-band fluxes. Column (8): Mid-infrared data available. 5.9, 6.75, 7.7, 9.63, 12.5, 14.3 and 16 fim indicate ISO data, and 10 fJ,m 
indicates ground-based iV-band measurements. Column (9): references for the mid-infrared fluxes. 

References in Columns (4). (6), (7) and (9): 1. Ho & Peng (2001). 2. This work. 3. Alonso-Herrero ct al. (2001). 4. Quillcn et al. (2001) 5. Clavcl et al. (2000). 6. Perez 
Garcia & Rodriguez Espinosa (2001). 7. Rigopoulou ct al. (1999). 8. Alonso-Hcrrcro ct al. (1998). 9. Maiolino ct al. (1995). 10. Rickc (1978). 11. Edclson ct al. (1987). 
12. Rudy, LcVan, 8z. Rodriguez- Espinosa (1982). 13. Lc Floc'h, private communication (2002). 14. McAlary, McLaren, & Crabtree (1979). 15. Bcndo ct al. (2002). 16. 
Neugebaucr ct al. (1979). 17. Scovillc ct al. (1983). 



Table 2 

Summary of high resolution optical, near-infrared and mid-infrared observations for additional 

Seyfert galaxies in the CfA sample. 



Galaxy 


Type 


WFPC1/2 


Ref 


NICMOS 


Ref 


KL band 


Mid-infrared 


Ref 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 3031 


Syl.5 


F547M 


1 


F160W 


2 


9 


10 /^m 


10 


NGC 3079 


Sy2.0 






F160W 


2 


9 


10, 16 /im 


6, 7 


NGC 3185 


Sy2.0: 












10, 12 fim 


11, 4 


NGC 4258 


Syl.9 


F547M 


1 


F110W, F160W, F222M 


3 


3 


10.5, 12.5, 17.9 fj,m 


3 


NGC 4501 


Sy2.0 










9 


10 Aim 


6 


NGC 4579 


Syl.9/L1.9 


F547M 


1 






4 


6.75, 9.63A*m 


5 


NGC 5194 


Sy2.0 






F110W, F160W, F222M 


2, 4 


8 


9.55, 12.5 Aim 


4 


NGC 7479 


Syl.9 


F569W 


1 


F160W 


2 


4 


6.75, 10, 14.3 Aim 


4, 9 


NGC 7743 


Sy2.0 






F160W 


2 


4 





Note. — See Table 1 for an explanation. Column (2): the Seyfcrt types arc from Ho ct al. (1997). Column (8): the mid-IR fluxes of NGC 4258 are 
ground-based. 

References in Columns (4), (6), (7) and (9): 1. Ho & Peng (2001). 2. Quillen et al. (2001) 3. Chary et al. (2000). 4. This work. 5. Clavel et al. (2000). 
6. Devereux (1987). 7. Perez Garcia & Rodriguez Espinosa (2001). 8. Ellis, Gondhalekar, & Efstathiou (1982). 9. Willncr et al. (1985). 10. Rieke & 
Lebofsky (1978). 11. Maiolino ct al. (1995). 



14 



Table 3 

Distribution of type Is and type 2s in different 

SAMPLES. 



Sample 


No. galaxies 


Type 1 


Type 2 


(1) 


(2) 


(3) 


(4) 


All b/a 


Extended CfA 


58 


43% 


57% 


RSA 


60 


28% 


72% 


Palomar* 


44 


27% 


73% 


12 fim. 


55 


45% 


55% 


b/a > 0.6 


Extended CfA 


47 


43% 


57% 


RSA 


41 


37% 


63% 


Palomar* 


32 


28% 


72% 


12 fim 


39 


56% 


44% 



Note. — Distribution of type 1 (1 — 1.5) and type 2 
(1.8 — 2) Scyfcrts based on their optical classification. 
* Only galaxies with certain Seyfert classification arc in- 
cluded. 
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Table 4 

IRTF K'L-banb Aperture Photometry and L-band unresolved contribution within 

A 3"-DIAMETER APERTURE. 



Galaxy 

(1) 


FWHM 

(2) 


Aperture 
(3) 


K' 
(4) 


K' -L 
(5) 


/unresolved 

fill 
J 3" 

(6) 


L-band Morphology 
(7) 


Mkn 334 


0.78 


1.5 


11.64 


1.53 


89% 


Compact 






3 


11.21 


1.44 










6 


10.97 


1.32 










9 


10.83 








Mkn 335 


0.72 


1.5 


10.88 


1.76 


94% 


Compact 






3 


10.60 


1.76 










6 


10.50 


1.72 










9 


10.46 








UM 146/UGC 1395 


0.84 


1.5 


13.19 


1.21 


65% 


Compact 






3 


12.63 


1.10 










6 


12.25 


0.82 










9 


12.00 








Mkn 590 


0.86 


1.5 


11.96 


2.25 


76% 


Compact 






3 


11.38 


1.95 










6 


11.00 


1.68 












10.80 








Mkn 993 




6 


11.40 


0.36 




Blind offsets 


NGC 1144 


1.11 


1.5 


13.19 


1.14 


40% 


Extended 






3 


12.22 


0.84 










6 


11.51 


0.60 










9 


11.05 


0.52 






UGC 6100 




6 


12.23 


0.23 




Blind offsets 


NGC 3227 


0.81 


1.5 


10.87 


1.37 


71% 


Compact 






3 


10.24 


1.15 










6 


9.82 


0.91 










9 


9.62 


0.76 






NGC 3362 




6 


12.80 


0.40 




Blind offsets 


NGC 3786/Mkn 744 


1.08 


6 


10.74 


1.07 


56% 


Tracking problems 


NGC 3982 




6 


12.67 


1.17 




Blind Offsets 


NGC 4051 


0.78 


1.5 


10.54 


1.59 


73% 


Compact 






3 


10.08 


1.44 










6 


9.80 


1.30 










9 


9.67 


1.21 






NGC 4235 


1.05 


1.5 


12.26 


1.03 


58% 


Extended 






3 


11.44 


0.76 










6 


10.77 


0.46 










9 


10.37 


0.24 






NGC 4253/Mkn 766 


0.69 


1.5 


11.37 


1.87 


90% 


Compact 






3 


10.96 


1.75 










6 


10.70 


1.60 










9 


10.55 


1.50 






NGC 4395 




6 


14.00 


1.29 




Blind Offsets 


NGC 4579 


1.02 


1.5 


11.28 


0.49 


45% 


Extended 






3 


10.32 


0.33 










6 


9.57 


0.22 










9 


9.17 


0.15 






NGC 5033 


0.98 


1.5 


11.59 


0.66 


35% 


Extended 






3 


10.59 


0.42 










6 


9.77 


0.26 










9 


9.31 


0.19 







16 



Table 4 — Continued 



Galaxy 


FWHM 


Aperture 


K 


K 


— L 


/unresolved 
f 3 " 

(6) 


L-band Morphology 


(1) 


(2) 


(3) 


(4) 


(■' 


5) 


(7) 


NGC 5256/Mkn 266SW 


— 


1.5 


13.21 





.93 


16% 


Extended 




3 


12.30 





.98 










6 


11.71 





.94 










9 


11.44 










NGC 5256/Mkn 266NE* 


— 


1.5 


13.51 





.58 




Extended 




3 


12.67 





.65 










6 


12.08 





.67 










9 


11.77 










NGC 5273 


1.6? 


1.5 


12.98 





.64 


38% 


Extended 






3 


11.92 





.48 










6 


11.17 





.27 










9 


10.83 










NGC 5283/Mkn 270 




6 


11.26 





.22 




Blind Offsets 


Mkn 461 




6 


L = 11.97 








Blind Offsets 


NGC 5347 




1.5 


12.80 


2 


.09 


65% 


Compact 






3 


12.09 


1 


.80 










6 


11.62 


1 


.51 










9 


11.34 


1 


.30 






NGC 5548 


- ■ 


1.5 


11.26 


1 


.93 


65% 


Compact 






3 


10.77 


1 


.82 










6 


10.46 


1 


.66 










9 


10.29 


1 


.53 






NGC 5674 


0.83 


1.5 


12.69 


1 


.41 


57% 


Compact 






3 


12.09 


1 


.30 










6 


11.66 


1 


.12 










9 


11.43 


1 


.03 






Mkn 841 


0.90 


1.5 


12.88 


2 


.10 


57% 


Compact 






3 


12.34 


2 


.04 










6 


12.10 


1 


.97 










9 


11.98 


1 


.92 






NGC 5929 


1.50 


1.5 


13.07 





.43 


16% 


Extended 






3 


12.13 





.27 










6 


11.45 





.04 






NGC 5930* 


1.95 


1.5 


12.37 





.63 




Extended 






3 


11.41 





.63 










6 


10.86 





.43 










9 


10.57 





.28 






NGC 5940 




1.5 


13.95 








Not detected 






3 


13.11 














6 


12.52 














9 


12.18 










Mkn 1513 


0.77 


1.5 


11.23 


1 


.68 


75% 


Compact 






3 


10.89 


1 


.67 










6 


10.78 


1 


.66 






UGC 12138 


0.86 


1.5 


12.39 


1 


.36 


54% 


Compact 






3 


11.81 


1 


.21 










6 


11.46 


1 


.07 










9 


11.29 





.98 
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Table 4 — Continued 



Galaxy 

a) 


r Wxllvl 

(2) 


Aperture 
(3) 


is* 

(4) 


IS* T 

(5) 


/unresolved 

(6) 


L-band Morphology 
(7) 


IN GO 7469 


0.75 


1.5 


11.17 


1.96 


7 l/o 


Compact nucleus + 






3 


10.41 


1.64 




Extended 






6 


9.90 


1.42 










9 


9.73 


1.34 






NGC 7479 


0.75 


1.5 


12.78 


1.59 


61% 


Compact nucleus + 






3 


11.95 


1.22 




Extended 






6 


11.19 


0.83 










9 


10.74 


0.64 






NGC 7603/Mkn 530 


0.69 


1.5 


10.85 


1.59 


100% 


Compact 






3 


10.45 


1.48 










6 


10.21 


1.37 










9 


10.07 


1.27 






NGC 7682 




6 


11.92 


0.05 




Blind Offsets 


NGC 7743 


1.00 


1.5 


11.93 


0.16 


< 45% 


Extended 






3 


11.22 


0.19 










6 


10.67 


0.20 










9 


10.38 


0.25 







Note. — Column (1): Galaxy. Column (2): Size (FWHM) of the nucleus as measured from the L-band image. 
Column (3): Diameter in arcsec of the circular apertures used for the photometry. Column (4): iiT'-band magnitude. 
Column (5): measured K' — L color. Column (6): Fraction of unresolved emission within a 3"-diameter apeture in 
the L-band. Column (7): Morphology of the L-band emission. 
* NGC 5930 and Mkn 266NE arc not part of the CfA sample of Seyfert galaxies. 



Table 5 

Galaxies with circumnuclear radio and mid-infrared emission. 



Galaxy 


6 cm i 13 
(2) 


10 M m / ~ 


Extended L 


Extended ISO 


(1) 


(3) 


(4) 


(5) 


Mkn 334 


0.12 


1.18 


no 


no 


NGC 1144 


0.07 




yes 


yes (at 6.75 fim) 


NGC 3982 


0.03 


0.34 




yes 


NGC 4051 


0.07 


0.77 


no 


no 


NGC 4388 


0.08 


1.90 


no 


yes + bright point source 


NGC 5033 


0.15 


0.19 


yes 


yes 


NGC 5256/Mkn 266SW 


0.10 


0.33 


yes 


yes + bright point source 


Mkn 841 


0.14 


1.28 


no 


no 


NGC 5929 


0.40 


0.41 


yes 


yes 


NGC 5940 




0.45 




yes 



Note. — Column (1): Galaxy. Column (2): Ratio of the nuclear (0.3") to extended (15") emission at 6 cm. 
Column (3): Ratio of ground-based to ISO fluxes at 10 (im. Column (4) and Column (5) indicate whether there is 
evidence for significant extended emission from the IRTF L-band images and ISO 9.63 fim images, respectively. 
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Table 6 

Nuclear (unresolved) optical and infrared fluxes for Seyfert 1.8 - 2 Galaxies in the extended CfA 

sample, and fitted spectral indices. 



Galaxy 


/opt 


/i i 


/l 6 


h i 


h 5 




fio 


fie 


CKopt — IR 


CKIR 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(9) 


(9) 


(10) 


(11) 


Mkn 334 


< 3.8 


2.33 


6.83 


13.8 


31.7 


131. 


148. 


< 390. 


1.80 ± 0.11 


1.80 ± 0.11 


Mkn 573 


< 0.097 


0.15 


0.54 


3.2 


18.8 




167.* 


< 550. 


2.98 ± 0.20 


2.91 ± 0.23 


UM 146/UGC 1395 


0.090 


0.35 


0.82 


3.0 


4.4 


15. 


15.* 


< 250. 


1.98 ± 0.18 


1.76 ± 0.18 


NGC 1068 




9.8 


98. 


445. 


3080. 


13600. 


19600. 


< 46900. 


3.02 ± 0.31 


3.02 ± 0.31 


NGC 1144 


< 0.020 


< 0.11 


0.08 




3.4 


< 120. 


< 158.* 


< 460. 


3.3 


3.3 


NGC 3079 






< 0.1 




< 61 




92.* 


< 2070 


3.7 


3.7 


NGC 3185 














20.* 








NGC 3362 


< 0.031 


0.04 


0.05 




< 2.2 




12.* 


< 210. 


2.75 ± 0.20 


2.75 ± 0.20 


U V. ■ D1UU 


*c u.uoo 


/ n 17 

<, U.l I 


<. u. io 






42 


A fi 
4:0. 


<. lOU. 




^ A 


NGC 3786/Mkn 744 


2.2 


1.75 


3.25 


< 11.9 


17.7 




59.* 


< 190. 


1.22 ± 0.19 


1.65 ± 0.09 


NGC 3982 


0.009 


< 0.27 


0.34 




< 7.2 


< 62. 


19.* 


< 430. 


2.77 ± 0.19 


2.4 


NGC 4258 


< 0.019 


< 0.5 


0.9 


4.0 


20. 




100.* 




2.85 ± 0.30 


2.39 ± 0.21 


NGC 4388 




0.06 


0.71 




39.9 


265. 


267. 


< 1550. 


3.51 ± 0.42 


3.51 ± 0.42 


NGC 4501 












< 35 


17.* 








NGC 4579 


0.202 






< 14 


12.5 


104. 


113. 




2.30 ± 0.13 


2.35 ± 0.46 


NGC 5033 


0.89 


0.61 


3.22 


< 8.0 


8.6 


< 179. 


24.* 


< 1340. 


1.41 ± 0.20 


1.67 ± 0.22 


NGC 5194 




< 0.14 


0.19 


0.38 


< 34. 




110. 




3.78 ± 0.11 


3.78 ± 0.11 


UGC 8621/1335+39 












49. 


49. 


< 130. 






NGC 5252 


0.070 


0.20 


0.70 


1.0 


11.1 


15. 


27. 


< 45. 


2.03 ± 0.20 


1.94 ± 0.26 


NGC 5256/Mkn 266SW 




< 0.08 


< 0.15 




2.4 


< 125. 


30.* 


< 330. 


3.0 


3.0 


NGC 5283/Mkn 270 


< 0.013 


< 0.43 


< 0.30 




< 11.0 


< 11. 


17.* 


< 60. 


2.7 


2.7 


NGC 5273 


0.27 


0.68 


1.67 




3.0 


23. 


25. 


< 90. 


1.63 ± 0.11 


1.66 ± 0.15 


Mkn 461 




< 0.45 


< 0.46 




< 4.7 


< 12. 


12.* 




1.8 


1.8 


NGC 5347 




0.22 


0.97 


4.0 


14.3 




208.* 




3.06 ±0.18 


3.06 ± 0.18 


Mkn 471 


0.15 


0.18 


0.32 




3.0 




14.* 




1.82 ± 0.18 


2.08 ± 0.14 


NGC 5674 


0.11 


0.90 


2.62 


4.2 


7.5 


38. 


43. 




2.08 ± 0.16 


1.78 ± 0.11 


NGC 5695/Mkn 686 


< 0.039 


< 0.17 


< 0.25 








< 10.* 








NGC 5929 


< 0.023 


< 0.17 


< 0.25 




1.3 


< 41. 


24.* 


< 550. 


2.9 


2.9 


UGC 12138 


< 2.4 


1.86 


2.59 


5.3 


9.0 


40. 


54. 


< 250. 


1.64 ± 0.08 


1.64 ± 0.08 


NGC 7479 


< 0.006 




0.24 


< 3. 


9.0 


150. 


263.* 




3.85 ± 0.24 


3.58 ± 0.30 


NGC 7674/Mkn 533 




1.25 


5.0 


12.3 


53. 


259. 


344. 


< 860. 


2.36 ± 0.18 


2.36 ± 0.18 


NGC 7682 


< 0.04 


< 0.10 


< 0.13 




< 4.8 


10. 


12.* 


< 290. 


3.0 


3.1 


NGC 7743 






< 0.5 


< 8. 















NOTE. — This table contains all the optically classified Seyfert 1.8 — 2 galaxies in the 'extended' CfA sample. The measured unresolved flux densities 
in Columns (2) through (9) are listed in mjy. Columns (10) and (11) arc the fitted spectral indices (j v cx v~ a ) in the optical-infrared (a pt-iR, spectral 
range ~ 0.4 — 16 fim) and infrared (otr, spectral range 1 — 16 pm), respectively. If no optical measurement is available, we give the same value for 
ttopt-iR and oir . When the errors arc not quoted for the fitted spectral indices it means that the formal errors are infinity (those cases with only two 
or fewer detections). 

Column (9) * indicates small aperture ground-based Af-band measurement (see discussion in Section 4.3). The ground-based TV-band fluxes (derived to 
give monochromatic fluxes of Raylcigh-Jeans spectra) need to be corrected to v f v — constant by multiplying them by a factor of 1.22 (see also Edclson 
ct al. 1987, and Section 5). 

Additional nuclear fluxes at 4.8 pm (ground-based) — Mkn 573: 41.3 mjy; NGC 1068: 8245 mjy ; NGC 5252: 15.7 mjy; NGC 7674/Mkn 533: 108.6 mjy 
(Alonso-Hcrrero ct al. 2001). Additional nuclear fluxes at 12 pm (ISO) — NGC 3185: 20 m Jy; NGC 5194: 400 mjy; NGC 5695/Mkn 686: < 55 mjy. 
Additional nuclear fluxes at 12 pm (ground-based): NGC 4258: 165 mjy. Additional nuclear fluxes at 14.3 pm (/SO) — NGC 7479: 520 mjy. Additional 
nuclear fluxes at 17.9 pm (ground-based) — NGC 4258: 300 mjy. 
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Table 7 

Nuclear (unresolved) optical and infrared fluxes for Seyfert 1 — 1.5 Galaxies in the extended CfA 

SAMPLE. 



Galaxy 


/opt 


fi l 


fl 6 


f 2 i 


f'i 5 


fe 75 


ho 


/l6 


^opt — 1R 




(!) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7)° 


(9) 


(9) 


(10) 


(ii) 


IVlKIl ooO 


4 






32 8 


78 2 


100. 


om 
zUo. 




i 9i _i_ n i s 

1 .Ol Hi U. lo 


n -i- n i c; 


UGC 524/0048+29 














41.* 








Mkn 993 


0.039 






< 3.1 


< ii. 




18.* 


< 130. 


2.2 


1.6 


1 Zw 1 


2.20 








110. 




310.* 


< 640. 


1.64 ± 0.14 


1.07 ± 0.10 


Mkn 590 


0.80 






9.3 


38.0 


144. 


198. 


< 460. 


2.08 ± 0.14 


1.81 ± 0.17 


IVlKIl lz4o 














16. 


< 210. 






NGC 3031 


1.39 




13.4 


< 202. 


< 161. 




86.* 




1.42 ± 0.15 


1.1 


NGC 3227 


2.70 




10.6 


22.6 


46.7 


294. 


382. 


< 1220. 


1.80 ± 0.10 


1.97 ± 0.11 


NGC 3516 


2.20 




18.1 




117. 


264. 


369. 


< 770. 


1.71 ± 0.11 


1.51 ± 0.16 


NGC 4051 


1.60 






38.2 


73.5 


265. 


411. 


< 1440. 


1.91 ± 0.08 


1.62 ± 0.07 


W C C /1 1 r: 1 
IN uC 41D1 


^1 A 
01.4 


oy.u 


i ni fi 

lUO.O 


1 77 ^ 


325 




i /inn * 


A i on 


i 9fi -I- n ns 

1 .ZD HZ U.Uo 


i A a -I- n m 
l . / ± z ± ill U . Uo 


NGC 4235 


0.44 




3.69 


< 6.6 


9.0 




34.* 


< 170. 


1.51 ± 0.09 


1.30 ± 0.04 


NGC 4253/Mkn 766 


3.90 




20.0 


20.3 


58.0 


177. 


274. 


< 730. 


1.65 ± 0.07 


1.56 ± 0.08 


Mkn 205 


2.10 








13. 


30. 


58.* 


< 190. 


1.34 ± 0.06 


1.53 ± 0.10 


NGC 4395 


0.17 




0.85 


< 1.1 


< 2.4 




12.* 




1.29 ± 0.10 


1.5 


Mkn 231 


16.3 




80.1 




326. 




1420.* 


< 3960. 


1.73 ± 0.07 


1.59 ± 0.06 


Mkn 279 


4.2 








32. 


102. 


146. 


< 290. 


1.14 ± 0.06 


1.38 ± 0.11 


NGC 5548 


4.3 




15.0 


31.6 


46.2 


177. 


268. 


< 440. 


1.52 ± 0.09 


1.43 ± 0.10 


Mkn 817 


2.2 








55. 


142. 


259. 


< 490. 


1.78 ± 0.11 


1.43 ± 0.05 


Mkn 841 


1.8 






4.1 


12.0 


74. 


133. 


< 340. 


1.75 ± 0.17 


2.18 ± 0.12 


NGC 5940 


0.46 




0.93 


< 1.9 




< 46. 


26* 




1.54 ± 0.18 


1.8 


Mkn 1513 








19.7 


44.6 


85. 


140.* 




1.08 ± 0.12 


1.08 ± 0.12 


NGC 6104 


0.12 




0.13 








12.* 


< 120. 


1.76 ± 0.35 


2.5 


NGC 7469 


6.4 


16.2 


39.0 


67.8 


90. 


433. 


600. 


< 1290. 


1.67 ± 0.09 


1.57 ± 0.09 


NGC 7603/Mkn 530 


0.53 






21. 


74.0 


138. 


176.* 




2.29 ± 0.32 


1.35 ± 0.19 



Note. — This table contains all the optically classified Seyfert 1 — 1.5 galaxies in the 'extended' CfA sample. The measured unresolved fiux densities in 
Columns (2) through (9) arc listed in mjy. Columns (10) and (11) arc the fitted spectral indices (f u oc in the optical-infrared (a op t — ir, spectral 

range ~ 0.4 — 16 fim) and infrared (cki;r, spectral range 1 — 16 (im), respectively. If no optical measurement is available, wc give the same value for 
a opt-iR an d ckir. When the errors are not quoted for the fitted spectral indices it means that the formal errors arc infinity (those cases with only two 
or fewer detections). 

Column (9) * indicates small aperture ground-based TV-band measurement (see discussion in Section 4.3). The ground-based TV-band fluxes (derived to 
give monochromatic fluxes of Raylcigh-Jcans spectra) need to be corrected to v f u — constant by multiplying them by a factor of 1.22 (see also Edclson 
ct al. 1987, and Section 5). 

Additional fluxes at 5.5pm (ISO) — I Zw 1: 173 mjy, NGC 4151: 882 mjy, Mkn 231: 606 mjy. Additional fluxes at 7.7 fim (ISO) — I Zw 1: 221 mjy, 
NGC 4151: 1090 mjy, Mkn 231: 919 mjy. Additional fluxes at 12 pm (ISO) — NGC 4395: 7.6 mjy. Additional fluxes at 14.3 pm (ISO) — Mkn 1513: 
160 mjy. 



Table 8 

Parameters of Efstathiou & Rowan-Robinson 
(1995) and Efstathiou et al. (1995) torus models. 



Model 


e c 


TUV 


fr/r 2 


h/r 2 


T 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Torus 1 


45° 


1000 


0.05 


0.5 


1000 


Torus 2 


45° 


500 


0.05 


0.5 


1000 


Torus 3 


45° 


250 


0.05 


0.5 


1000 


Torus 4 


30° 


1200 


0.01 


0.1 


900 


Torus+Cone 


30° 


1200 


0.01 


0.1 


900 



Note. — Column (2): O c is the half opening angle of the 
torus. Column (3): t\jv is the equatorial UV optical depth 
in magnitudes. Column (4): r\jri is ratio of the inner 
and outer torus radii. Column (5): h/r 2 is the ratio of the 
height scale and the outer radius. Column (6): T (in K) is 
the dust sublimation temperature. 




Fig. 1. — Contours plots of selected galaxies with extended emission in the L-band (right panels). We show the if'-band contour plots for 
comparison (left panels). The orientation of the images is south up, east to the left. 
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Fig. 2. — Le/i panel: Non-stellar SED of NGC 1068 in units of ergcm~ 2 s _1 . The ISO spectrum of the AGN (nuclear) component of 
NGC 1068 is from Le Floc'h et al. (2001). The additional mid-infrared data points at 12.6, 17, 22.5, 24.5, and 33.5 fim are ground-based small 
aperture observations from Rieke & Low (1975a); data points at 11.2 and 20.5 fim are from Alloin et al. (2000). Right panel: Non-stellar SED 
of NGC 4151 in units of erg cm -2 s — . The additional mid-infrared data points at 10 and 18 fim (unresolved emission) are from Radomski ct 
al. (2002), at 21 ^tm is from Rieke & Low (1972) and at 33.5 /im from Rieke & Low (1975b). 
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Fig. 3a. — Observed non-stellar SEDs (y /„) of optically classified type 1 (1 — 1.5) CfA Seyfcrts in units of erg cm 2 s 1 . 
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Fig. 3b. — Continued. 
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Fig. 4a. — Observed non-stellar SEDs (y /„) of optically classified intermediate type (1.8 — 1.9) CfA Seyferts in units of ergcm~ 2 s~ 
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Fig. 4b. — Continued. 
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Fig. 5a. — Observed non-stellar SEDs (y f v ) of optically classified type 2 CfA Seyferts in units of erg cm 2 s 1 . 



27 









1 


1 




1 






1 1 11 






1 1 




O 
* — i 




Mkn461 






NGC5347 


A ~ 




NGC5695 






NGC5929 ^_ 





cv 




u 


J, i A 




A 

A 

A 

A 






I 4. 

I 


4- - 

I 




4- 

A 

A _ 

I 


































H — 1 1 1 1 II h- 


1 1 1 1 II h 




H — 1 1 1 1 II 1 — 


+4+] h 




H — 1 1 1 1 II 1 — 


++4+] h 




1 10 


"S 


o 

* — 1 

cv 




NGC7674 

A 

A 

A 


A A ^ _ 

A 




NGC7682 


i ~ 

A 




NGC7743 

i 

X 








GO 


' — 1 
















I 





































. — 1 































1 


10 




1 


10 




1 


10 







\ (yttm) \ (/2m) A (;ii,m) 



Fig. 5b. — Continued. 
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Fig. 6. — Distributions of the fitted spectral indices (/„ oc v~ a ) in the optical-infrared wavelengths (a pt-lR f° r the spectral range 
~ 0.4 — 16 fim, left panels) and infrared wavelengths (am for the spectral range 1 — 16 (im, right panels). The very top panels are the 51 
Scyfcrts in the 'extended' CfA sample with measured spectral indices. The panels below are the same sample broken up in different Seyfert 
types as derived from optical spectroscopy. 
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Fig. 7. — Effects of foreground extinction (Ay = 0.4, land 5 mag) on the observed SED of an AGN, using Rieke & Lebofsky (1985) 
extinction law. The AGN SEDs arc outputs of the Torus+Cone model from Efstathiou & Rowan-Robinson 1995 (see Section 6 and Tabic 8 
for details) for three different viewing angles. 8 V = 0° is the polar view of the AGN, 8 V = 30° is an intermediate view (corresponding to the 
half opening angle of the torus, C ), and 9 V = 90° is the equatorial view (through the torus). 
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Fig. 8. — Fitted optical-infrared (left panel) and infrared (right panel) spectral indices (see Section 5) as a function of the inclination of the 
host galaxy, where b/a = is an edge-on galaxy, and b/a = 1 is a face-on galaxy. The arrows at the top of the plots indicate the values of 
b/a for the remaining 7 galaxies with no measured spectral indices — dotted lines for the type 2s (green in the on-line version) and solid lines 
for the type Is (red in the on-line version). The dashed line indicates a host galaxy axial ratio of b/a = 0.6 or an inclination of i > 53° (see 
discussion in text). The galaxies are separated between Seyfert Is (1 — 1.5) shown as open circles (red in the on-line version) and Scyfert 2s 
(1.8 — 2) shown as filled star symbols (green in the on-line version), using the classification from optical spectroscopy. 
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Fig. 9. — Predicted infrared spectral indices (/„ oc v~ a iR) in the 1 — 16 fim spectral range as a function of the cosine of the viewing 
angle to the torus (8 V = 0° polar view of the AGN, 8 V = 90° equatorial view of the AGN) using Efstathiou & Rowan-Robinson (1995) and 
Efstathiou et al. (1995) torus models. Left panel: results for models with a torus half opening angle C = 45°, and two different equatorial 
UV optical depths: ttjv = 500 mag (filled triangles), ttjv = 250 mag (filled stars), Models Torus 2 and Torus 3 respectively (see Table 8). 
Right panel: results for a torus model with a half opening angle C = 30° , and an equatorial UV optical depth of ttjv = 1200 mag with the 
cone component (Model Torus+Cone in Table 8, open circles) and without the cone component (Model Torus 4 in Table 8, filled circles). 
In both diagrams the shaded areas represent the observed ranges of fitted «tr for Seyfert 1 — 1.5s (darkest shaded area), Seyferts 1.8 — 1.9s 
(intermediate shaded area), and Seyfert 2s (lightest shaded area) in the CfA sample. 



